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AN INTEGRATED THEORY FOR PREDICTING THE HYDRO- 
THERMOMECHANICAL RESPONSE OF ADVANCED 
COMPOSITE STRUCTURAL COMPONENTS 
by C. C. Chamis,* R. F, Lark,* and J. H. Sinclair* 

ABSTRACT 

An integrated theory is developed for predicting the hydrothermo- 
mechanical vHDTM) response of fiber composite components. This 
integrated theory is based on a combined theoretical and expc *1 mental 
Investigation. In addition to predicting the HDTM response of components, 
the theory is structured to assess the combined hydrothermal effects on 
the mechanical properties of unidirectional composites loaded along the 
material axis and off-axis, and those of angleplied laminates. The 
theory develoiied predicts values which are in good agreement with 
measured data at the miciomechanlcs, macromechanics, laminate 
analysis and structural analysis levels. 

Key Words; Com|X)sites, theory, cxix?riment, comjiosite mechanics, 
laminate theory, structural analysis, stress analysis, warpage, thermal 
stresses, hydro stresses, finite element analysis 


Aerospace Engineer, Natioi al Aeronautics and Sixice Administration, 
Lewis Research Center, Cleveland, Ohio 44135. 
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NOMENCLATURE 


A 

a 

b 

C 

D 

E 

F 

m 

G 

HDTM 

K 

k 

M 


m 




w 


x,y,z 

1,2,3 

a 

/3 


(Property and direction are ftlven by subscrlptB) 
inplane stiffness array (eq. (32)) 
dimension 
dimension 

coupled stiffness array (eq. (32)), heat capacity 
array of bending stiffness (eq. (32)) 
elastic modulus 
function 

mechanical property relation (eq. (4)) 
thermal property relation (eq. (5)) 
shear modulus 

abbreviation for hydrothermomechanical 

heat conductivity 

volume ratio 

moment 

moisture 

transformation matrix for strain and stress respectively 

ply uniaxial strength 

glass transition temperature 

resin softening temperature = T 

reference temperature (273 K) 

out-of- plane displacement 

structural axes reference coordinates 

material axes reference coordinates 

thermal expansion coefficient 

moisture coefficient 
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^sub 



K 

V 

p 

a 


f I 


{ } 

I 1= f 
Subscripts 


void magnification factor 

correlation coefficients, taken as unity if not defined 
strain 

reference plane inplane strain (eq. (32)) 
ply orientation angle taken jwsitlve for transformations 
from the old to the new axis, 
local curvature (eq. (32)) 

Poisson's ratio 

density 

stress 

strain magnification factor 

square array 

transpose 

inverse 

column vectt)r 

row vector 


C 

c 

d 

f 

K 

H 

f 

m 

r 


S 


T 

t 


c(mipression 

composite property 

dry property 

fiber property 

glass transition 

hydro (moisture) property 

ply (unidirectional) com|Xisite proi^erty 

moisture i roperty 

resin proix?rty 

shear property 

tension 

thermal pn.'perty 
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V void 

w wet property 

x,y,z property or variable referred to structural axes 

o reference temperature (T^ = 273 K) dry property 

1,2,3 property or va/iable referred to material axes 

INTRODUCTION 

'ihe influence of hydrothermal environments on advanced composite 
mechanical properties has been extensivelv investigated by the composites 
community over the past several years (refs. 1 to 5). The various in- 
vestigations thus far have generated significant fundamental knowledge 
for: (1) identifying the variables and mechanisms which bring about 
this influence. (2) establishing key relationships between these variables, 
and (3) providing the basis for developing an integrated theory to pre- 
dict the hydrothermomechanical (HDTM) resjxinse of advanced composite 
structural comixments. The objective oi this investigation is the develop- 
ment and description of an integrated theory for predicting the HDTM 
response of advanced comjxisite structural components. 

The approach jxirsued in developing this integrated theory con- 
sisted of the following: (1) review of work to date in order to select re- 
lationships between key hydrothermal variables (moisture and tempera- 
ture) of resins and in order to identify and/or develop simple approx- 
mate relationships for relating resin physical and mechanical prop- 
erties to the key hydrothermal variables, (2) performance of limited ex- 
|X?riments to verify the conclusions arrived at fron reviewing the wor’; 
in (1), (3) incorporation of the simple relationships (identified in the liter- 
ature or developed herein) as approp;.*iate into comjiosite micro- 
mechanic'’, comiJosite macromechanics, laminate theory, and compr)Site 
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structural analysis to complete the development of the desired inteKrated 
theory for the HDTM response, and (5) application of the intef^rated theory to 
selected examples and comparison with experimental c!ata. Some of the 
RoverninR equations are given in matrix notation in the text for con- 
venience. Their expanded form is given in the Appendix. 

EXPERIMENTAL INVESTIGATION 

The experimental program was conducted to assess the effects of 
water absorption on the weights and volumes of an epoxy matrix resin 
and fiber^epoxy resin composites as well as to provide data to compare 
with comiX)site micromechanic predictions. 

Materials Investigated 

Tbe materials evaluated in this investigation consisted of neat 
PR-288 epoxy matrix resin and an AS-tyije graphite fiber PR-288 resin 
unidirectional composite. 

Materials Fabrication 

The neat PR- 288 ejx^xy matrix resin was supplied by the 3 M 
Company in plate form. The fiber/resin laminates were made at the 
Lewis Research Center using commercial AS-tyjje graphite PR-28R 
(Ai^^'E) resin prepreg. Individual plies of prepreg were stacked in a 
metal mold to form 10 ply and 25 ply unidirectional laminates. A 
thermocouple was inserted into the laminate at one <^nd. The cold nx)ld 
was placed in a hydraulic laminating press previously heated to 450 K 
(350° F). The press platens were adjusted to apply contact pressure of 
0. 10 MPa (15 psi) on the mold. After the laminate temperature reached 
311 K (100° F), contact pressure was maintained for another 3 minutes. 

The platen pressure was then increased gradually over a 2 minute period 
to 2. 1 MPa (300 psi). The laminate was tlien cured at this pressure for 
a period of 2 hours. After completion of curing, the laminate was re- 
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moved from the preBS and the mold in a hot condition and allowed to cool 
at room temperature conditions. 

Specimen Pieparation and Testintt 

The neat epoxy resin and fiber resin composite specimens were cut 
into test specimens using a diamond cutting wheel. The dimensions and 
designations of the specimens are given below. 

Type and Dimensions Designation 

Neat Resin 

0. 64 cm (0. 25-in. ) cube Thick Resin 

1. 27 by 1. 27 by 0.32 cm Thin Resin 

(0. 50 by 0.50 by 0. 125-in. ) 

Fiber Resin Composite 

10 ply - 1. 27 by 1. 27 by 0. 13 cm Thin Composite 
(0. 50 by 0. 50 by 0.05-in. ) 

25 ply - 1 . 27 by 1 . 27 by 0. 25 cm Thick Composite 
(0.50 by 0.50 by 0. 10-in. ) 

Ten test specimens of each tyi)e were prepared and characterized for 
weight and physical dimensions in the as-cut condition. 

All of the specimens were placed in a 339 K (150° F) oven to remove 
absorbc'd water and were allowed to remain in the oven until weights and 
dimensions were stabilized. Tlic specimens were then subjected to room 
temperature [294 to 297 K (70° to 75° F) and approximately 50 percent 
relative humidity] and boiling water environments. The weights and di- 
mensions of the specimens were periodically determined. 

Results and Discussion 

Figure 1 shows the effects of a boiling water moisture environment 
on the percentage changes in specimen volumes and weights for the thin 
epoxy resin specimens. Tliese changes increased nonlincarly with time. 
Similar percentage changes for the thick epoxy resin specimens were 
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erratic and are not shown in figure 1. This erratic behavior is attributed 
to possible nonuniform moisture distribution through the specimen thick- 
ness. 

Fixture 2 shows the effects of a iKjiling water moisture environment 
on the percentage chant; es in specimen volumes and weitthts for the graphite 
fiber epoxy resin composite specimens. TTie percenlat;e chanttes for both 
the thick and thin composites increase nonlinearity with ti- The percent 
chanttes in volume and weight for the thin specimen were essentially 
stabilized at the termination of the tests. However those of the thick 
composite were still changing gradually. 

Figure 3 shows the relationships between the percentage changes in 
the volumes and weights of the specimens suujected to room temperature 
and lK)illng water environments, respectively. Tbe slopes of the lines 
drawn through these data are equal to the coefficients listed in table I. 

Tlie percentage dimensional changes versus weight changes (moisture ex- 
pansion coefficients) are summarized in table II. 

The points to be noted from figures 1 and 2 are that the absorbed 
moisture appears to have reached equilibrium in the thin resin and thin 
composite specimens, but not in the thick composite specimens. The 
moisture content in the resin at equilibrium is about 7.8 percent which is 
comparable to that reported in the literature (ref. 1). As can be seen 
both weight and volume changes are the same. The corresponding changes 
for the thin composite are about 3 percent. Tins is equivalent to about a 
7 percent change in the amount of resin in the composite which is about 
40 percent by volume. 

The points to be noted from tables I and II are that the coefficients from 
the room temperature exjjosure specimens are 1 for the resin specimens 
and 0.3 for the composite specimens. The coefficients from the boiling 
water exi)osure are about 1.1 for all the specimens. This indicatss that 
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the volume ex|)aruiion of the resin in the composite is the same as in the 
next resin and also that the in situ moisture density is about the same as 
that of water. 

THKORY, COMPAHIiiONS, AND DISCUSSIONS 
In this part of the paper are described the flow diagram of the inte- 
grated theory for hydrothermomechanical (HDTM) response, the concepts 
and assumptions for developing this theory, the key governing equations 
for composite micromechanics, composite macromechanics, and laminate 
theory. Finally, the incorporation of all of these into a computer code 
and the coupling of this code with a structural analysis procedure for pre- 
dicting the HDTM response of advanced composite structural components 
are described. Comparisons are made with available data. 

Flow Diagram of Integrated TTieory for 
Hydrothermomechanical (HDTM) Response 
The flow diagram of the integrated theory for HDTM response of 
composite components is illustrated schematically in figure 4. Tlie key 
elements in this diagram are blocked with single lines while the desired 
result is blocked with double lines. As can be seen in the diagram the 
key elerients are composite micromechanics, composite macromechanics, 
laminate analysis, and structural analysis. Tlie desired result is the 
HDTM response of composite components. Each key element receives 
input from the element preceding it and supplies input to the element 
following it as well as additional information. In a design procedure the 
HDTM response must be compared to the appropriate design criteria in 
order to assess the adequacy of the composite component. 

Concepts and Assumptions 

The tw'o concepts underlying this integrated theory arc based on: 

(1) the observation that "the moisture affects the resin HDTM properties 
only, " (refs, 1 to 5) and (2) the hypothesis that the moisture effects on 
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the resin HDTM properties can be expressed as a function of the use 
temperature and the reduced (;lass transition temperature of the wet 
resin The reduced ({lass *.*ansition temperature is 

related to both ({lass transition temperature of the dry resin 
and the amount of moisture in the resin (m^.) as described in reference 6. 
This hypothesis Implies iso|)arametric functional relationship between 
the various resin properties and may be expressed in symbolic form as 
follows: 

resin mechanical properties T^^. | (1) 


resin thermal properties (o, K. » Fjliof , K, T^^ | (2) 


where E,v, and S denote modulus, Poisson’s ratio and fracture stress, 
respectively; and where a, K and C denote thermal expansion coefficient, 
heat conductivity and heat caiwcity, respectively. Moisture diffusivity and 
electrical conductivity are expected to have similar functional relationships. 
The ({lass transition temperature of the wet resin T^^. is {{Iven in 
symbolic form as follows: 


^({wr “ ^3 ^^ndr’ '^r^ 


(3) 


where, as mentioned previously, T^dr’ *^r ({lass transition 

temperature of the dry resin and m^ the moisture content in the resin. 

(Sec ref. 6 for the functional form of Fg). 

The followlnR assumptions were made in developin({ this inteprated 
theory: (1) the moisture diffusion in the composite is independent of the 
stress state, (2) the temperature distribution in the composite is independent 
of the stress state. (3) all the assumptions that are common to composite 
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linear micromechanics and macrumechanics, and linear laminate theory; 

(4) the moisture effects in the composite are manifested through the resin 
as described by micromechanics; and (5) the hydrothermomechanical 
response of the composite structural component may be determined by 
solving the appropriate moisture diffusion, heat transfer, nd structural 
problems as will be indicated later. 

Relationships for Wet Resin Properties 
The relationships among the wet resin mechanical properties and the 
room temperature dry properties and temperature (eq, (1)) were deter- 
rrdned in this investigation usit)g the following procedure: (1) retention 
ratios (wet property at test temperature divided by room temperature 
dry property) of wet resin and unidirectional composite properties from 
the literature (refs. 1, 4, 7, and 8) were plotted versus temperature, 

(2) a sim|ile algebraic expression wa*. sought to approximate the retention 
ratio and thereby establish the functional relationship Fj in equation (1). 
The motivation for seeking a simple algebraic exjiression is that this 
can readily be incorporated into available composite micro mechanics 
equations. Tlie procedure used to establish the relationship is illustrated 
in figure 5. As can be seen the simple algebraic relationship 

IT - T ~ ^ ^ 

wet resin mechanical ^n^ert^ at test tempera ture _ ^ | s 

dry resin mechanical property at room temjierature ^s " ^o 

(4) 

is a good approximation of the data in the range 200 K (-100^ F) < T < T 

O 

Tlie notation in equation (4) is as follows: denotes mechanical prop- 
erty retention ratio, Tg equals T is the use temperature and T^ 

equals 273 K (0“ C). Tlie temperatures are expressed in K. The moist- 
ure effect is incorpoiated through T^ as defined in equation (3). Tlie 

I I r ’ I — 
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corresponding relationship for the wet resin thermal properties Is assumed 
to be the inverse of the equation (4) or 


wet r esin therma l property at te st temperature ^ 
dry resin thermal property at room temperature ^ 


fT - 

T 

B 

0 

- 

T 

L * 



1/2 


(5) 


Equation (5) Is a good approximation for the literature dry resin data 
(refs. 9 and 10) In the range 200 K (-100° F) < T < as is shown in 
figure 6. It is anticipated that equation (5) will fit wet resin thermal prop- 
erties equally as well in the same range. If It does not,, another relation- 
ship may he found and the procedure described below can be used to in- 
corporate the moisture effects thermal properties. Summarizing then; 




ro 


m 


( 6 ) 


(a,K.C)^^»(a,K,0^„ (7) 

where the subscrl|it ro refers to dry resin pro|)erties evaluated at 
room temperature; and and are given by equations (4) and (5), 
respectively. 

Composite Micromerhanics 

The governing micromechanics equations which incorporate the 
moisture effects in the unidirectional composite (ply) response were 
derived using the procedures of references 11, 12, and 13. Since the 
number of these equations can be very large (refs. 14 and 15) herein 
we present those equations which are either used to predict results or 
are used to illustrate important phenomena. 

Volume an d density relatio nships . - The ply density (pf^) and ply 
moisture (m ^) are given by; 

"fw = "(d * <'>fd * * '‘v* >"r 


( 8 » 
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IK, ♦!(,) 

(9) 

where; 


^fd * ‘^f^f + **r ^r 

(10) 


(11) 

1 a kj + k j. + k^ 

(12) 


and where the notation in equations (8) to (12) is as follows; p denotes 
density; k, actual volume ratio; m, moisture; moisture expansion 
coefficient; the subscripts f, r, v, and f denote fiber, resin, void, and 
ply pro|X'rty, respectively; and the subsc.*<ots d and w dimotc dry 
and wet con JiMons, respectively. The presence of voids affects all the 
proper^*..., given by e<iuations (8) to (11). Note that equation (11) can 
Le used to predict the theoretical moisture expansion coefficient 
of the resin by setting = 0. For example, typical resins have 
specific gravities of about 1.2. The theoretical volumetric exiansion 
coefficient of those resins is 3*3^ which is equal to the si>eciflc gravity of the 
resins or 1. 2 in this example and is in good agreement with the measured 
coefficient of 1 . 1 of the boiling exposure specimens in table I. TTie value of 
of about 0.4 is also in good agreement with literature data (ref. 1). 
Equation (9) predicts a moisture content of alxmt 2.6 |x>rcent for an 
AS PR-288 comix)Site com|)ar‘xi to about 2.9- percent for the thin com- 
|X)6ite and 2. 3-|X?rcent for the thick comjx)8ite from figure 2. These 
com|)arisons show’ that the values predicted by equations (9) to (11) 
are in good agreement with measured data. They also show thnt the 
resin volume expansion is conserved in the composite as the composite 
exi^ands at these moisture contents. This means tliat the resin is essen- 
tially incomjiressible and will probably behave like a viscoelastic or 
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vlscoplaatlr medium. A direct conclusion from this condition is that 
microstress relaxation occurs in the composite. Another conclusion is 
that the voids in the In situ resin may collapse at the high (about 
1 -percent and greater) moisture content. 

Ply elastic constants . - The governing micro mechanics equations 
for predicting the ply elastic properties longitudinal modulus 
transveise modulus (£^ 22 ^* Poisson's ratios (>^^ 12 ) intra- 
laminar shear modulus respectively, are; 

(13) 

E 

£^22* F= 

^mEro/Ef22) 


'fl2 " ^fl2 


+ *'r 'Vo 


(15) 


fl2 


a” G 
m ro 


1 -iT?;(i - 


’"mGro/Gfi2) 


(16) 


where 


' ro 

"2(> ♦ >'ro> 


( 1 .) 


and where is given by equation (4); the subscript f denotes 

fiber property; the subscripts 1 and 2 denote fiber directions with 1 ta’;en 
alo.Tg the fiber and 2 transverse to the fiber direction. Note that the major 
Poisson's ratio is independent of moisture as may be verified from equa- 
tion (17). 

P ly e;q)ansivi ties. - The thermal expansion coefficients («fjj 
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and o f22^ moisture expansion coefficients or swelling coefficients 

and respectively, by the following micromechanics 

equations: 


"ni 



^^ro ^m*^r •*'°fn*^f^fll 
^ m^ro’^r ®^fll‘'fll 


(18) 


®f22 * “ro’^t^' ■ *^f *'ro®fll/®'m^ ■*■ ®f22^^ 

^fll ' ^r*^r '^m®'rc/®^fll 


(19) 

( 20 ) 


^C22 ” 






^ro 


|/lCf E j22 ■*■ (' ~ t ^ 


( 21 ) 


m^ro 


or for Incompressible resin with /3^jj « 0. 



3li^kj.p 


fd 


'<3^Wp^ 


( 22 ) 


where is given in equation (4) and in equation (5). At small 

moisture (less than 0.3 percent) contents ft (02 probably closer to 
the value predicted by equation (21) than equation (22) which was derived 
assuming conservation of volume of the in situ resin. 

Ply uniaxial strengths. - The uniaxial ply strengths, tranverse 
tension (Sj22t^ compression (Sj22C^ intralaminar (in- plane) 

shear (Sfj2S^ given, respectively, by the following micromechanics 
equations: 
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Sj22t 


®f22C 


^fl2S 



f22 





^fl2 


(23) 


(24) 


(25) 


where denotes void magnification factor and <p^ strain magnification 
factor, S denotes fracture stress; the subscripts T, C, and S denote 
tension, compression and shear, respectively. The coefficient ^22T' 
^22c» coefficients if needed, otherwise they 

are assumed to be unity. The void and strain magnification factors are 
given by: 


/3^ = 1/1 - (4k/7Tkj.)j 


1/2 


^/i22 ' 


1 kf Ef22'^®' 


f ^ f 22^ ^f22 

1 


(26) 

(27) 


^ni2 ;z 

1 - vkf 


(28) 


The uniaxial ply longitudinal tensile streng^th ^S^nj) and compressive 
(Sjiic) are givem by: 

SfiiT =Sj.j, +kj. 


(29) 
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®mc * |[®fC ^^CS®fl2S *®rCo^l 

where S^j and denote fiber longitudinal tensile and compressive 
strength, respectively; the coefficients and are correlation 
coefficients If needed, otherwise taken equal to unity; is given by 
equation (4) and S^|2S equation (25). It is important that 
the various correlation coefficients in equations (23), (24), (25), 

(29), and (3) be evaluated from dry room temperature data and 

should remain invariant with both temperature and moisture changes. 
These coefficients as >.ell as ^CS in equation (30), are evaluated 
using the procedure described in reference 12. Though a specific set 
of micromechanics equations were used herein, the functional relations 
and should work equally well with any other set. 

Comparisons with measured data . - The unidirectional (ply) prop- 
erties predicted by micromechanics for an AS/3501-5 composite, using 
the fiber and resin properties given in table m are tabulated together 
with literature data (ref. 8) in table IV. The properties arc for one 
moisture content (1.8 jjercent) and two temperatures (room and 366 K 
(200° F)). By comparing corresponding measured and predicted prop- 
erties it is seen that the agreement is about 10-percent, which is con- 
sidered very good, for most of the properties except at room 

temperature, Sj.j2s at both temperatures, and E^^22 ^fl2 ^t 366 K 

(200° F)). One explanation may be the normal difficulties encountered 
in measuring these properties. There is no measured wet data for the 
thermal expansion coefficients. 


I 


! 
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The predicted value for from equation (20) with ^i^fc0.41, 

kj,*0,4 and 0205 is about 0. 04 compared with 0. 05 

for the thin composite and with 0.02 for the thick composite, table II. 

This difference may be due to possible irreversible danuq^e induced in 
the resin during cooling down. The predicted value for ^^22 
equation (22) with k^ = 0, and Pf^j = 1.58 is about 0.79. This is in 
good agreement (within 2-perccnt) with the averages of the thickness 
and width moisture expansion coefficients for the thin specimen, 0.75, and 
for the thick specimen, 0,80 (table II). 

The significant result from these compaiisons is that the micro- 
mechanics equations predict ply MDTM properties which are in good 
agreement with measured data at different hydrothermal environments. 
Tliis agreement leads to the conclusion that the micromechanics equations 
given herein may used to assess the composite response in different 
hydrothermal environments for preliminary design. Another significant 
conclusion is that the moisture resistance of various resins including 
new ones may be assessed by measuring the room temperature prop- 
erties given in table III and tliose properties required to determine T . 

s 

Co mposi te Mac ro mechanics 

The governing macromechanics equations pertinent to this discussion 
arc Uie ply HDTM relatitinships and the ply combined-stress failure 
criteria. TIu' ply MDTM relationships in matrix form are given by 

u f |(7f I + ^Tf {nr^} + } (31) 

where j } denotes a 3 x 1 vector and (1 a 3 ^ 3 array, their cxjiandetl 
form is indicated in the appendix. AT^ and m^ denote the ply tern- 
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perature and moisture respectively. Equation (31) may be solved for 
Oj if needed. Equation (31) is written about the i^y material axes (1,2,3) 
with 1 along the fiber 2 transverse to it, and 3 through the thickness. 
Equation (31) may be transformed to any other axes by well known trans- 
formations (refs. 10, 14, 16). Note the similarity between ATj(Oj) 
and m^ This similarity offers the following two features: (1) these 

terms arc conputationally interchangeable, and (2) they may be combined 
on an element by element basis, 'Hiese two features make it possible to 
use available coni|)uter codes which include only the term to pre- 

dict the HDTM response of angleplied laminates. Reference 14 provides 
such a computer code and an integrated theory for the structural re- 
sponse of angleplied laminates subjected to thermal and mechanical 
loads. This code can also be applied to the i^itegrated theory for the 
HDTM response of angleplied laminates by incorporating: (1) the two 
above features and (2) the moisture effects through the micromechanics 
equations described previously. 

The ply combined-stress failure criterion used in this investigation 
is that develoixjd in-house (LeRC combined stress failure c/lterlon which 
is described in ref. 12) and is available in the computer code (ref. 14). 

Laminate Analysis 

Governing e quatio ns. - The laminate analvsis jart of the integrated 
theory for predicting the HDTM resjxinse of comix>site components is 
available in the comjxiter code (ref. 14), The governing laminate analysis 
equations include those for predicting the HDTM resjJonse of the laminate 
and those for predicting the ply HDTM strains and stresses. The 
equations for the laminate HDTM response are: 
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where N, M, e, k are 3 v i vectors (see appendix for expanded form) 
denoting, respectively: inplane force, bending moment, reference plane 
strains, and local curvatures; A^, C^, anci D^, are 3 v 3 arrays 
denoting; axial, coupled, and bending stiffness, respectively. These arrays 
are determined from ply properties using laminate analysis procedures 
(refs. 8, 14, and 16). The subscript c denotes composite (laminate) 
property or variable, and the subscripts m and t denote moisture 
and thermal properties. Note that and are the thermal 

forces and moments and and arc the hydro forces and 

moments. Equation (32) is an important result of laminate theory and 
is generally used to predict the following: 

1. HDTM laminate relationships required in structural analysis. 

2. Laminate forces (displacements) when the displacements (forces) 
are known from structural analysis 

3. Ply HDTM stresses 

4. Residual and hydro stresses 

5. Laminate warjxige from hydrothermal stresses. 

The ply HDTM strain and stresses are, respectively, given bv: 

=|R,cKI<-co' I'cl' 

|o,l = |E,K|n,^| |f,| - m, 1(3,1 - iT, (a,)) (34) 

where is a 3 x 3 strain transformation matrix which is a function 
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of 9 measured from the laminate structural axes to the ply material 
axis (see appendix for expanded form), and is the normal distance 
from the reference i^ane to the center of the ply. The ply hydrothermal 
strains and stresses for symmetric laminates with uniform temperature 
and moisture conditions are, respectively, given by: 

- (o^)) (35) 

\ - {(if!) + AT(fR^JIa^) - (a^))) (36) 

where the subscript HT denotes hydrothe~mal property and the other 
symbols have been previously defined. As can be observed from equa- 
tions (35) and (36) combination of moisture (m) and temperature (AT) 
may exist which will result in zero hydrothernul strains and stresses. 
As can also be observed hydrf)thcrmal environments when both moisture 
and teni|Jeraturc increase or decrease induce the largest hydrothermal 
strains or stresses. 

Application laminate elastic properties. - The laminate analysis 
just described was used to analyze some of the laminates tested in ref- 
erence 8. Tlie results obtained for elastic properties are summarized 
in table V. The comparisons show that the predicted values are in good 
agrwment with measured data (within 10 percent) except for the shear 
moduli of laminates 1 and II where the predicted values are 18 and 
15 percent higher, respectively. A possible reason for this large dif- 
ference is inaccuracy of the experimental values arising from the dif- 
ficulty in measuring shear properties. 

Application^ hydrothcrmome£hamcal stresses. - The predicted ply 
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HOTM stresses for one of these laminates are summarized in bar-chart 
form in fit^ure 7, both individually and combined. The uniaxial strenKth.s 
are also shown. Tlie following; points are worthy of note in figure 7: 

1. The hydro and residual stresses have opposite si(mB with the 
hydro stress dominating (about 3-times as high for some stresses). 

2. The hydro and mechanical stresses have the same sign for the 
longitudinal and transverse stresses in the 45^ plies and the trans- 
verse stress in the 0° plies. 

3. The mechanical load is the dominant source of the shear stress in 
the 45^ plies. 

4. The transverse residual st»*ess magnitude is alx)ut equal to the 
uniaxial transverse strength in the 0° and 45° plies. It is ]x)ssible that 
this stress caused traiujply cracking in these plies prior to moisture con- 
ditioning of the laminate. 

5. The mechanical load shear stress in the 45° plies exceeds the 
corresponding uniaxial strength 

6. The hydrothermal stress contribution to the longitudinal stress 
in the 0° plies is negligible compared to mechanical stress. 

7. It appears that fracture was initiated by shear failures of the 
45° plies and completed when the 0° plies failed in longitudinal tension. 
One ini|X)rtant conclusion from the above discussion is that the ply 
hydrothermal stresses in angleplied laminates are substantial and 
need to be considered in both analysis and design. 

App lic ation - laminate fractuj^e. - Predicted laminate fracture 
stresses for combined HDTM load conditions arc comjxired with 
measured data (ref. 8) using a bar-chart in figure 8. Note in this 
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fi(^re both measured data rani^es and predicted rani^es are shown by 
horizontal lines In the bars. The lower bound of the predicted fracture 
stress is based on first ply failure (usually Induced by Iransverse or 
shear stress) while the upper bound is based on failure of the remaining 
major load carrying plies (not necessarily fiber fractures). The pre- 
dicted fracture stresses are based on room-temperature-dry uniaxial ply 
strengths nnd room-temperature-wet ply hydrothermal elastic properties. 
This was done for two reasons: (1) the lower bound on laminate fracture 
stress predicted using room-temperature-wet ply strengths was con- 
sidered to be unrealistically low and, (2) to compe'^sate for possible 
in situ ply transverse and intralaminar shear strengths which are 
suspected to be higher than the uniaxial values tested under the same 
HDTM environment. 

Examining the bars in figure 8, it is seen that the measured and 
predicted ranges overlap in 13 of the 15 cases. The two cases where 
they do not overlap are laminate 1 under longitudinal tension and 
under longitudinal compression. In both of these cases the predicted 
values are lower than the measured data by about the same amount. 

One r(>ason the predicted longitudinal compression fracture stress 

is lower may be the in situ enhancement of the ply longitudinal compressive 

strength. 

Tlie important observation from the data in figure 8 is that dry ply 
strength data should l)e used to predict laminate fracture stresse.s in 
order to ol)tain reasonable agree*ment with the measured data. Another 
one is that the predicted low'er bound on fracture stress appears to be 
on the conservative side. Tlie results of this stress analysis show that 
the LeRC ply comiiined-stress strength criterion can be used for pre- 
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dii'tinit laminate fracture under combined HDTM loads. One significant 
conclusion from the previous discussion is that the in situ ply strengths 
may be hi(;her than those determined under uniaxial tests. 

Application - hydrostresses. - Theoretical ply transverse hydro- 
stresses in (I'/O/ - ^/0)g laminates of different composite systems are 
shown in fi^re 9 as a function of 0 . These hydrostresses corresjxind 
to 2. 5 percent moisture and were predicted using equation (36) with 
AT = 0 and m = 0.025, As can be seen from the curves in figure 9 the 
transverse hydrostresses are compressive, increase with increasing 
M and leveling at alx)ut 0 = 60°, and are sensitive to composite system, 
ranging from about 0.034 GPa (5 ksi) for Kevlar/epoxy (KEV/E) to alwut 
0.19 GPa (28 ksi) for boron epoxy (B/E). The corresponding intra- 
laminar shear stresses arc shown in figure 10, The intralaminar shear 
stress increases very rapidly with 0 peaking at about - 45°. It is 
sensitive to composite system ranging from about 0.014 GPa (2 ksi) for 
KEVE to 0.083 GPa (12 ksi) for B E. 

Comparable plots in |±" |g laminates are shown in figure 11 for 
transverse stress and in figure 12 for intralaminar shear. In these 
laminates the ply transverse hydrostress increases very rapidly, with 
0 peaking at 45°. and is symmetric about 45°. It is strongly influenced 
by the composite system ranging from ,*bout 0.034 GPa (5 ksi) for KEV/E 
to alx)ut 0. 19 GPa (28 ksi) for B E, Note that these ranges are about 
the same as those for the \<^/ 0 / - f^/0|g except that the variation with ft 
is different. The corresponding intralaminar shear stresses are shown 
in figure 12. This stress is antisymmetric with respect to 45°. has two 
peaks at about 0 30° and 0 60°, and depends on the composite sys- 

tem, ranging from alx)ut 0,014 GPa (2 ksi) for KEV/t; to alxjut 0.069 GPa 
(10 ksi) for B/E or S -glass /epoxy (S-G/E). 
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The major conclusion to be drawn from the above hydrostress results 
is that the transverse stresses are compressive in an IncreaeiuRly uni- 
form moisture environment and as such they may enhanc e the ply in situ 
transverse strength. 

The ply hydro stresses induced in a B/E unidirectional laminate 
during moisture absorption are shown in figure 13 and during de- 
sorption In figure 14. The moisture profiles used in these predictions 
were obtained from reference 1. ITte curves show that large trans- 
verse stresses can be induced. These stresses are compressive 
in the outer plies and tensile in the inner plies during moisture ab- 
sorption. During desorption the transverse stress is tensile in the 
outer and center plies and compressive in the interniodiate plies. 

Tlie longitudinal stress is negligible for both cases. Note that curves 
for both stresses have the same but reflected shape of the moisture 
curve. The magnitude of the transverse tensile stress is about 0.090 
GPa (13 ksi) for the absorption case and is about 0.048 GPa (7 ksi) for 
the desorption case compared to the corresponding strength of about 
0.055 GPa (8 ksi). One conclusion fn)m these results is that rapid 
moisture absorption or desorption may induce transply cracks. I- or a 
more extensive discussion on hydro stresses see reference 17. 

Structural Analysis 

The structural analysis ixu-t of the integrated theory described 
herein consists of the coupling of the composite mechanics computer 
code (ref. 14) with a general purpose finite element structural analysis 
program, for example. NASTRAN (ref. iG). Referring to figure 4 the 
integrated theory for the HDTM response is completed when the 
laminate analysis supplies the material properties (eq. (32)) required 
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In NASTRAN to generate the element atlffne^e, and when NASTRAN pre- 
dirtes the HDTM atructural reeix^nse displacement vector (c^q and 
oq. (32)) required by laminate theory to perform the laminate 
stress analysis. 

A simple example of HDTM structure response of a composite 
component Is the free warpa^e of an angleplied laminate subjected to 
nonuniform moisture dlstributicn through the thickness. One way 
to assess this warpage is to fix the laminate at one corner and deter- 
mine the lateral displacement at the diagonally- opposite corner. This 
structural response may also be predicted by the following close 
form equation fur a rectangular flat laminate 

" ' \ ' 2 (37) 

where w is the corner displacement, a and b the side dimension, and 

Kp are the local curvatures which are determined from equation (32). 

Corner displacements obtained for a square laminate a = b = 25.4 cm 

(a = b = 10 in. ) for two laminate configurations |(i30)«| and (30/0/ - 

^ s 

30/0 Ig from an AS/E composite and subjected to three different non- 
uniform through-the- thickness moisture p.ofiles are summarited in fig- 
ure 15. As can be seen in this figure the corner displacement due to 
warpage induced by moisture can be substantial. Note that corner dis- 
placements greater than 2.54 cm (1-in. ) are probably beyond the limits 
of linear structural analysis. It is interesting to note that the linear and 
paral)olit moisture profiles induce approximately the same waipage. 

The variation of the warpage corner displacement as a function of ply 
angle is shown in figure 16 for laminates with hyperbolic moisture pro- 
files. Note that the warpage is sensitive to laminate configuration and 
is sensitive to certain moisture profiles. One conclusion from this dis- 
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cusBlon is to take advantage of these sensitivities to assess stress relaxation 
due to viscoelastic (viscoplastic) behavior of moist resins. 

CONCLUSIONS 

Tlie major results and conclusions of an investigation to provide an 
Integrated theory for predicting the hydrother mo mechanical (IIDTM) 
response of advanced composite components are as follows: 

1. An integrated theory was developed to predict the HDTM resj)onse 
of adx’anced composite components. This theory Incorporates the hydro- 
thermal effects In composites through micromechanics, macromechanics 
and laminate analysis all of which are lurts of a com|x>slte mechanics 
comixiter code. The comjjuter code is coupled with a major structural 
analysis finite element cominitcr program to complet3 the integrated 
theory. 

2. A simple relationship was established which relates the mechanical 
properties of the wet resin to its room temjxjrature dry proix?rties and 

to its wet glass-transition temperature. The corresiX)nding thermal 
properties relationship is the reciprocal of the mechanical properties 
relationship. 

3. The resin conserves its volume in a hydro environment. The 
density of the in situ moisture is about the same as water. Voids in the 
in situ resin may collapse at high moisture contents (greater than 1- 
percent). 

4. Composite micromechanics predicts hydro expansion and ply 
mechanical pn)perties (elastic and strength) which are gencrc.lly in good 
agreement with measured data under different hydrothermal environments. 

5. Laminate theory predicts laminate wet properties which arc 
within 10-percent of measured daUi. 

6. Laminate theory in conjunction with the LeFtC failure criterion 
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predicts laminate fracture ctre^s ranges which overlap with tnnse of 
the measured data. The predicted lower bound on fracture stress 
tends to be below the measured data scatter in general. 

7. The hydro stresses in the plies of angleplled laminates reach 
magnitudes in the transverse direction and in intralaminar r*hear which 
are comixirable to corresponding ply strengths. 

8. The ply hydro transverse stress is either tensile or compression 
depending o:i whether moisture is absorbed or desorbed. The profiles 
of ply stresses through-the-thickness of the laminate are mirror images 
of the corresponding moisture profile. The ply hydro stresses depend 
on moisture content and profile, laminate configuration and composite 
system. KFV/E angleplied laminates have the lowest ply hydro stress 
and n^E or S-G E the highest. 

9. Nonuninform moisture through the laminate thickness induces 
severe laminate wari>age. This warjiage depends on lamii.ate configuration 
and can prcJuce corner displacements several times the laminate thickness. 

10. Hydro stress in angleplied laminates can be predicted using the 
equivalance feabire between hydro and thermal expansions. 

11. In sit«' ply strength may be enhanced by hydrothermal environments 
as compared with uniaxial data. The in situ resin may behave like a 
viscoelastic or viscoplastic moHium in hydrothermal environments. 

Moisture induced warpage in angleplied laminates may be an effective 
means to assess stress relaxation in a hydrofnermal environment. 
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APPENDIX 

(Vector and array expansion) 



I'fi-Kfii- ‘ 

f22’ ^112^ 


and similarly for all other vectors. 



1 1/E,„ 


0 

' ■‘'fl2^® f22 

1/E ^22 

0 


0 

0 

1/Gfj2 ’ 


^cll ^cl2 

^cl3 


fAcI = 

^c21 ^c22 

^c23 



^c31 ^c32 

^c33 



and similarly for fE^. |, [C^ ), and [D^ ]. 

{Np } - |Np I - Npyy, N^J^yl 

and Similarly for } , I } , and 1 } . 

cos^f) sin^P -Sin20 
“1 T 2 9 

f^fcjl ^ cos'^fl sin2P 


^sin20 --sin2P cos20 
2 2 
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TABLE I. - MOISTURE PrXPANSlON 


VOLUME COEFFICIENTS FOR 
EPOXY RESIN ANB FIBER RESIN 


COMmSITE SPECIMENS 


Specimen 

Coefficientj 

♦ 

Room temperature 

1 

1 

1 

Thin resin 

1.0 

Thick resin 

1.0 

Thin comixjsite 

.3 

Thick comjjositc 

.3 

Boiling water 


Thin resin 

1.0 

Thick resin 

1.1 

Thin comiX)Site 

1.2 

Thick comixisite 

1.1 


TAMLK II. - MOISTl'HK FXPANSION COKFUl'IKNTS 
FOK LINFAL DIMFNSIONS FOR FIMFH RFSIN 
COMPOSITF SPFCIMFNS 


Specimen 

Coefficients 


Spe<’imen directions 


LentJth 

Width 

Ttiickness 

Averaue thickness 
width 

Tliin composite 

0.05 

0.5 

1.0 

0.79 

Thick composite 

0.02 

0.9 

0. 7 

0.80 
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TABLE m. - rroER AND RESIN PROPERTIES AND COEFEICIENTS fSEB 


IN THE MICROMECHANICS EQI'ATIONS 


Property 

Symbol 

I'nita 

Value 



St (Enxllah) 


AS fcraphllr liber 




Lonffttudlnal modulua 

^Tll 

UPa (MSI) 

221 (32) 

Tranaverar modulua 

*’F22 

CPa (MSII 

13.8 (2.0) 

Shear modulua 

^'F12 

GPa (MSI) 

13.8 (2.0) 

Pnlvion'a ratio 

L iKitudInal thermal exp. mefl. 
^"BpaverBe thermal exp. coefl. 

‘FI 2 

•fii 

“f22 

"k‘’('’f‘’) 

o^-I(Of-I) 

0.2 

•WIO"* (-0. 56^10'*) 
l»I0‘* (5.6v|0*®) 

'-ngltudlnal tenalle atreoKth 

®FT 

GPa tKSI) 

2.7 (400) 

..onKitudInal rompreaalon 

Sfc 

GPa (KSI) 

2.7 (400) 

3501 -S Ream room temperature 
(rel. 8) dry 




Modulua 


GPa (MSI) 

4.1 (0.4) 

0.4 

5. 7'10**(32^10‘*’ 

Thermal exp. eoell. 
Mnlature exp. roell. 

•'ro 

®ro 

"K‘’("f ') 

0.4 


Tenalle alrenRth 

®rTo 

GPa IKSI) 

0.048 (7.0) 

Compreaalon atrenitth 

®rCo 

jPa (KSI) 

0.25 (36.3) 

Shear atrenKth 

*rao 

CPa 'KSI) 

0.048 (7.0) 

Coelfteienta 




equation (30) 

— 


12.1 
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TABLE IV. - COMPARIJ^ON OK MEASURED AND PREDICTED MEC HANICAL PROPERTIES OF ANGLEPLIED LAMINATES 


AS/3501-5 IN HYDROTHERMAL ENVmONMENTS 


Property 

Symbol 

Units 

SI 

(English) 

Properties with 1.8 percent moisture 

Room temperature 

200® 1 

• 

Measured 

Predicted 

• 

Measured 

Predicted 

Density 

P( 

kgm (lb in.^) 


1.6'10^ (0.57) 

— 

1.6 10^ (0.57) 

Lunttitudinal modulus 


GPa (MSI) 

131 (19.1) 

133 (19.4) 

135 (19.5) 

133 (19.3) 

Transverse modulus 

^f22 

GPa (MSI) 

8.6 (1. 25) 

7.25 (1.05) 

4.21 (0 61) 

5 18 <0 75) 

Shear modulus 

1 

^fl2 

GPa (MSI) 

3.73 (0.54) 

3.52 (0 51) 

1.73 (0 25) 

2 28 (0.33) 

Poisson’s ratio 


— 

0.30 

0.28 



0 28 

Temp, expansion coeff. 

•’ni 

K-1 (Op-1, 


-5.6'10*^ (-0.31 -lO*®) 



-4. 0*10'^ (-0.22-10*®) 

j Temp, expansion coefl. 

”(22 

K*1 (Op-1, 


2 3'10*® (0.13 10*^) 


2 60 10*^ (1.55*10'^) 

1 Loni:. compressive strength 

^fllC 

GPa (KSI) 

0.891 (130) 

1.11 (16: ) 

0.891 (130) 

0.319(1341 j 

1 Lomt. tensile strength 

®fllT 

GPa (KSI) 

1.65 (240) 

1.66 (242) 

1 f6 (242) 

1.65 (241) 1 

1 Transverse tensile strength 

®f22T 

GPa (KSI) 

0.034 (5.0) 

0.038 (5.5) 

0 020 (2.9) 

0 023 (3.3) i 

Transverse conpressive strength 

^I22C 

GPa (KSI) 

0.202 ( 29.5) 

0.181 (26.4) 

0.110 (16.0) 

0. 119 (17 3) 

In-plane shear strength 

®M2S 

GPa (KSI) 

0.045 (6.6) 

0.033 (4.8) 

0 035 (5. 1) 

0.023 (3 8) ^ 

Thickness 

t, 

cm (in. ) 

0.014 (0.0055) 

0.0140 (0.00552) 


0.0140(0.00552) ! 


1 





. ^ 


From reference 8. 
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TAHLK V. - CONtPAHISON OK MKASl Kri) ANU PRKDICTKt} 
KLASTIC PHOPKHTtt:S OK AN(il,K:PLtK.lJ IJkMINATKai 
AS 3501-5 WITH 1.8 PK;WK:NT MOISTIHE 


Utmlnate 

Lome. 

modulus 
OPa IMSIl 

Trane, 
modulus 
OPa (MSI) 

Shear 
modulus 
GPa (MSI) 

Major 

Poisson's 

ratio 

I'O -45, O *451 

* a 

‘Measured 
Predlrted 
■t difference 

43 (6.3) 
43 (6.3) 
0 

21.2 (3.08) 
23 (3.2) 

• 3.9 

22.2 (3 21) 

26.2 (3,80) 
♦ 18.4 

5.54 (0.803) 
5.39 (0 781) 
-2.70 

irOj i45 Oj 90 O^ 

‘Measured 
Predicted 
(i difference 

89.7 (13.0) 
89.7 (13.0) 
0 

29 (4.2) 
31 (4.5) 
♦ 7.1 

10 (1.5) 

11 (1.6) 
• 6.7 

2.24 (0.325; 
2.19 (0.3181 
-2.2 

IU|(0 i45 '90),| 

* a 

‘Measured 
Predicted 
X difference 

46. 1 (6.68) 
49.7 (7.20) 
♦ 7.8 

45.7 (6.62) 

49.7 (7.20) 
♦ 8.7 

16.1 (2.34) 
18.6 (2.70) 
♦ 15.4 

2.42 '0.350) 
2.30 (0.333) 
-4.8 


'Room temperature wet property, MeaauretJ lUU from reference 8. 



Tiq. 1 • fAiislure absorption characteristics of thin epory resin specimens in 
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Fig 3. - Comparison of weight and M>lume changes of resin and composite speci- 
mens In boiling water. 
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Ft? ' . - CODMirHon ot triMSurNl «nd prMlicM fracture stresses of anqie- 
r'<' (animates AS/3501-5iMtH 1 8% moisture. 
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Fiq. 13. - Predicted stresses induced by moisture in unidirectian4i 
boron/epoxy leminete exposed to itmosphere of lOOl relative 
humidity for 28 days 
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F 19 . 16 - PredicteO warpaq* induced by T. 5 percent moisture in AS(E 
laminates with one surface protected (Hyperbolic moisture profile, t 
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